Majorana fermions in density modulated p-wave superconducting wires 
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We study the p-wave superconducting wire with a periodically modulated chemical potential and 
show that the Majorana edge states are robust against the modulation of density. In some particular 
regimes, the Majorana fermions are greatly enhanced by the periodic modulation. The existence 
of Majorana edge fermions in the open chain can be characterized by a topological Z2 invariant of 
the bulk system, which can be applied to determine the phase boundary between the topologically 
trivial and nontrivial superconducting phases. We also demonstrate the existence of the zero-energy 
peak in the spectral function of the topological superconducting phase, which is only sensitive to 
the open boundary condition but robust against the disorder. 
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Introduction.- Searching for Majorana fermions (MFs) 
in condensed matter systems has attracted intensive 
studies in past years [1, 2], due to the fundamental inter- 
est in exploring the new type of particles fulfilling non- 
Abelian statistics and the potential applications for the 
topological quantum computing [1-6]. Among various 
proposals for realizing the emergent MFs, the quantum 
wires with p-wave pairings provide promising candidates 
for realizing emergent MFs at the ends of wires. As shown 
in the original work of Kitaev [3] , boundary MFs emerge 
in the one-dimensional (ID) p-wave superconductor if 
the system is in a topological phase. Moreover, quan- 
tum wires with a strong spin-orbit coupling, or topologi- 
cally insulating wires, subject to a Zeeman magnetic field 
and in proximity of a superconductor, are found to ex- 
hibit boundary MFs Particularly, experimental sig- 
natures of MFs in hybrid superconductor-semiconductor 
nanowires have been reported very recently [9], which 
stimulates the study of exploring MFs in ID systems 0- 
[l^. Schemes of realizing Majorana chains have also been 
proposed in cold atomic systems [17], carbon nanotubes 
[1^, superconducting (SC) circuits [l9| and quantum- 
dot-superconductor arrays [13]. 

As most of theoretical works focus on ideal homoge- 
neous wires with a uniform chemical potential, an im- 
portant problem is the stability of the MFs under the 
modulation of density and disorder. In this work, we ex- 
plore the nonuniform p-wave SC wires with a periodic 
modulation of the chemical potential, extending the Ki- 
taev's p-wave SC model [s*], and study the fate of the 
MFs under the density modulation and disorder. In gen- 
eral, a periodically modulated potential is of benefit to 
the formation of periodic density waves. If the modula- 
tion amplitude is large, one may expect that the SC phase 
could be destroyed. As the robustness of Majorana mode 
is protected by the SC gap, the MFs may be unstable in 
the presence of density modulation and disorder. On the 
other hand, recent work on the density modulated wires 
in the absence of SC order parameter indicates the exis- 



tence of topologically protected edge states [2l|, |22| . So 
far, it is unclear whether the Majorana edge states is 
enhanced or suppressed under the density modulation. 
The current work will focus on this problem and show 
the robustness of Majorana edge states against the den- 
sity modulation and disorder. Particularly, in some pa- 
rameter regimes, we find that MFs always exist for the 
arbitrarily strong modulation strength, which may shed 
light on the design of quantum architectures producing 
robust MFs. 

Model of density modulated p-wave superconductor s.- 
We consider a typical lattice model of the ID p-wave su- 
perconductor with modulated chemical potentials, which 
is described by 

H = ^[(-tcjci+i + Aqq+i + H.c.) - iiic\ci], (1) 

i 

where c\ (ci) is the creation (annihilation) operator of 
fermions at the i-th site, t the nearest-neighbor hopping 
amplitude, A the p-wave SC order parameter and chosen 
real, and the modulated chemical potential given by 



Mi 



Vcos(27rza + 6) 



(2) 



with V being the strength, a = p/q a rational number {p 
and q are co-prime integers), and S an arbitrary phase. 
The modulated chemical potential can be generated by 
the bichromatic optical lattice for cold atom systems 
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23[ or through the control of gates for quantum-dot 
arrays and quantum wires [23, 113 • 

In order to diagonalize the quadratic form Hamilto- 
nian, we resort to the Bogoliubove-de Gennes transfor- 
25 , 26| and define a set of new fermionic opera- 



mation 
tors, 

1 ^ 

Vi = 7, y][(^n,i + V^n,i)4 + (^n,i - V^n,i)Ci], (3) 



where L is the number of lattice sites and n = 1, • • • , L. 
For convenience, (j)ri,i and are chosen to be real due 
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to the reality of all the coefficients in the Hamiltonian ([T]) . 
In terms of the operators r]ri and r]}^ , the Hamiltonian ([T]) 
is diagonalized as = Yln=i ^nivliVn — where 
is the spectrum of the single quasi-particles. From the 
diagonalization condition, [r]n,H] = A^/^n, we can get 
the following coupled equations: 

An^n,i = (A-t)V^n,i+l -MiV^n,i - (A + t)V^n,i-l, 
An^n,i = -(A + t)^n,i+l -/^i(/>n,i + (A-t)(/)n,i-l, 

(4) 

with z = 1, ••• ,L. From Eqs. one can prove 

that if the equations have the solution of (^n,*, V^n,*) 
{i = 1, • • • , L) with a positive eigenvalue A^ > 0, 
(^n,i, —i^n,i) is also the solution with the eigenvalue 
-An, which implies r]n{^n) = vU-^n)- 

Boundary Majorana fermions.- The solution to Eqs. 
(|4|) is related to the boundary condition. The boundary 
Majorana states can only exist in the system with open 
boundary conditions (OBC). To seek such a state, we 
solve Eqs. (j4]) under OBC of V^n,L+i = V^n,o = ^n,L+i = 
0n,o = 0, and the Majorana edge states correspond to 
the zero mode solution of A^ = 0. For A^, = 0, 0^ and i/ji 
are decoupled and Eqs. (|4]) can be written in the transfer 
matrix form 

with Ai = I j . To understand the boundary 

Majorana mode, we rewrite the operators (|3j) as 

L 
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a=l/2 
v=i (3) 



a=l/3 
V=1 (d) 



1 ^ 



(6) 



where = c] -\- Ci and = i{ci — c]) are operators 
of two MFs, which fulfill the relations (7f )^ = 7f and 
{7f ,7^} = 2SijSa(3 with a and f3 taking A or B. Coeffi- 
cients of (j)ri,i and tjjn^i in Eq. (|3j) are just the amplitudes 
of Majorana operators and 7^^, respectively. If there 
exists a zero mode solution, according to Eqs. (|5]), we 
will get a decaying solution for one set of coefficients and 
a growing one for the other. An example is the spe- 
cial case with V = and A = t for which the zero 
mode solution is given by (^1, ^2, • • • , 0l) = (1, 0, • • • , 0) 
and (V^i,-- - = (0, • • • ,0,1). In this case, 

r]^{A = 0) = I [75^ + ^7l]' which means that the zero 
mode state is divided into two separated MFs located at 
the left and right ends, respectively. 

For the commensurate potential with a = p/q^ Ai is q- 
periodic. One can judge whether the zero mode solution 
exists by evaluating the two eigenvalues of the product 
matrix A = Yli^i^i (HI- If both of them are either 
smaller or greater than unity, there is a zero mode so- 
lution with MFs located at the ends, otherwise no zero 
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FIG. 1: The energy spectra of ID modulated p-wave super- 
conductors vs the phase factor, S, under OBC for a = | 
((a)-(c)) and | ((d)- (f)). The red dashed lines are the lowest 
excitation states and the blue dotted ones are the edge states 
in the excitation gaps. These figures are calculated with pa- 
rameters t = A = 1 and L — 101. 



mode solution is available. Solving Eqs. (j4]), we can 
get the whole excitation spectrum for the single quasi- 
particles. In FigdJ we show the spectra for cases of 
a = 1/2 and 1/3 under OPC. It is shown that the zero 
mode solutions exist in the whole phase parameter space 
when V is smaller than a critical value. When V ex- 
ceeds a critical value, the regimes for the existence of the 
zero mode solutions in the phase parameter space become 
narrow. However, our numerical results indicate that the 
zero mode Majorana solutions always exist at some spe- 
cific values of 5, no matter how strong the strength V is. 
For the systems with a = 1/2 and 1/3 shown in FigdJ 
the specific values of 5 are (^^+^)^ for a = 1/2 and 
(2m+i)7r ^ _ with m being an integer. At these 
specific values Sg^ the critical value of the phase transi- 
tion, Vc{A^Ss)^ which is a function of A and J, tends to 
infinity. When deviated from these points, there exists a 
finite critical value, Vc{A^S)^ above which there appears 
a topologically trivial phase without zero mode MFs. 

In order to see clearly the differences between the zero 
mode state and nonzero mode states, we display distri- 
butions of and i/ji for the zero mode solution and solu- 
tions with nonzero eigenvalues of quasi-particles located 
at the bottom and center of the continuous band (orange 
squares along the green cut in Fig. ffl^c)). As shown in 
Fig 121 foi" the system with S = Stt/S and V = 3, the 
amplitudes of the Majorana operators for the zero mode 
solution are located at the left and right ends, respec- 
tively. As (i/ji) decays very quickly away from the left 
(right) edge, there is no overlap for the Majorana modes 
of and . As a comparison, distributions of and 
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a=1/2 



FIG. 2: The amplitudes of 0^ (black solid) and ipn (red 
dashed), for excitation states labelled by orange squares along 
the green cut in Fig. [U^c). From bottom to top, (a), (b), and 
(c) represent the states of zero mode, at the bottom and at 
the center of the bulk band, respectively. They are calculated 
with parameters y = 3, t = A = l, and 

L = 101 under OBC. 



tpi for nonzero modes spread over the whole regime and 
the corresponding quasiparticle operator can not be split 
into two separated Majorana operators. 

Z2 topological invariant.- The presence or absence of 
zero mode MPs is determined by the Z2 topological class 
of the bulk superconductor. As no boundary zero mode 
solution is available for the system with periodic bound- 
ary conditions (PBC), we can choose a Z2 topological 
invariant (or 'Majorana number') [3] to characterize the 
topological nature of the bulk system. In order to define 
such a topological invariant, we shall consider the sys- 
tem with PBC. For the periodic system with a = p/q^ 
it is convenient to use the Fourier transform, Ci = Cs^i = 
Y^Xl/c Cs,/ce*^'^^ with i = s -\- {I — l)g, which transforms 
the Hamiltonian ([T]) into H = Z]/c[Z]sIi (~^<^I,/c^s+i.fc + 
Acs,kCs^i-k) - tcl f^ci^ke'^^ + Acq,feCi -fee"'^^ + H.c] - 

5^s=i /^s(4,fc^s,/c- |), where 5 = 1, • • • , is the index 
of inner sites in a supercell, / = 1, • • • ^L/q the index of 
the l-th supercell, and k the momentum defined in the 
reduced Brillouin zone of [0^27r/q). Then we define a set 
of new operators [12]: 



72s-l(^) = Cs,k + 



72s (A:) = {cs^k 



(7) 



with the anticommutation relations {7^1(^)5 7n (^01 = 
2Smnhk' and jl,{k) = jm{-k). 7m(0) and jm{7T/q) 
are just Majorana operators due to 7^(0) = 7m (0) and 
7m(^/^) 7m{-7T/q) jm{-7T/q)' In the basis of the 
new operators, we can transform the Hamiltonian into 
the form: 



H = j^^Bm,n{f^hm{-k)-fn{k), 



(8) 
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FIG. 3: Z2 topological invariant, M, vs the phase factor, S, 
for system with t = A = 1, a = ^ ((a)-(c)) and | ((d)-(f)) 
under PBC. 



with B2s-i,2s{k) = -B2s,2s-i{k) = -/is for s 
l,---,q, B2s-ias+2(k) = -B2s+2as-i(k) = A - t, 
B2s,2s+i{k) = -B2s^i,2s{k) = A + t for 5 = 1, • • • ,g - 1, 
Bi,2q{k) = -Bl^ ,{k) = -(A + t)e-*^^ and B2,2q-i{k) = 
-B^2,-ia{k) = -{A-t)e-^^^. 

The parameters Bm,n{k) form di 2q x 2q matrix B{k), 
and here only B(0) and B{7r/q) are skew-symmetric. Fol- 
lowing Kitaev [3|, ll2| , we can calculate the Z2 topological 
invariant defined as: 



M = sgn[Pf(5(0))]sgn[Pf(5(^/g))], 



(9) 



k rn,n 



where Pf(X) = EpSgn(P)Xp,P2 ' ' ' 

the Pfaffian of the skew-symmetic matrix X with P 
standing for a permutation of 2N elements of X and 
sgn(P) the corresponding sign of the permutation. Ac- 
cording to the definition, generally we have M = ±1 with 
M = 1 corresponding to a Z2-topologically trivial phase 
and M = — 1 to a Z2-topologically non-trivial phase. In 
Fig. 3, we display the Z2 topological invariant versus the 
phase factor 6 for systems with the same parameters as 
in Fig. [TJ Comparing Fig. 3 with Fig HI we see the exact 
correspondence between the presence (absence) of zero 
mode MPs and —1 (+1) value of the Z2 topological in- 
variant. Especially, M takes the value of (blue square 
dots in Fig.3(b)) at F = 2, which is just the critical 
point Vc of the phase transition from a Z2-topologically 
non-trivial phase to a topologically trivial phase for the 
corresponding S. In principle, we can always determine 
Vc{A^S) through the condition of M = and give the 
whole phase diagram. 

Phase diagram.- Without loss of generality, we choose 
A, F > in the following discussion. From the def- 
inition of the matrix B{k)^ we know that for a = 
1/2, Pf[P(0)] = -V^ cos^ S - and Pf[P(7r/2)] = 
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FIG. 4: Phase diagrams expanded hy V, A, and 6 for the cases 
of (a) a = 1/2 and (b) a = 1/3. The curved surfaces are the 
phase boundaries above which is Z2-topologically trivial while 
below which is Z2-topologically non-trivial. 

-V^cos^5 ^ 4A2. Here due to Pf[5(0)] < 0, we just 
need to set Pf[5(7r/2)] = to get the phase bound- 
ary condition, that is, V'^ cos^ 6 = 4A^, or equivalently 
V\cosS\ = 2A (Fig. m^a)). Likewise, for a = 1/3, 
Pf[5(0)] = cos3S - 2t{f + 3A2) and Pf[5(7r/3)] = 

— ^ cos3J + 2t{t^ + 3A^). The phase boundary condi- 
tion can be also got easily as V^ \ cos3(5| = 8t(t^ + 3A^) 
(Fig. m^b)). From Fig. |4l we see clearly again that there 
exist some specific points 5^, at which the system is al- 
ways Z2-topologically non-trivial for arbitrary V. On the 
other hand, the critical value, Fc(A, (5), increases with the 
increase of the SC pairing amplitude A. 

For the general case of a = p/g, we can in- 
fer the phase boundary condition by analyzing the 
expression of Pfaffian. The non-permutated term, 

Bi2Bs4 - - - B2q-l,2q, glVCS (" 1)"^ fls^l A^s ' the 

sparsity of the matrix 5, we see that there are only 
three kinds of permutations which contribute non-zero 
terms: 1. 52s-i,2s^2s+i,2s+2 ^2s-i,2s+2^2s,2s+i or 
^1, 2^2^-1, 2g ^i,2(?^2,2g-i, which makcs a replace- 
ment of /ig/ig+i or /ii/iq in the non-permutated term by 
(A^ -t^) or (A^ -t2)e-2^^^ respectively; 2. B2s-i,2s ^ 

^2s,2s+l {S = l,---,g-l) and B2q-l^2q ^ ^l,2g, 

which gives -{t + A)^e-*^^; 3. B2s-i,2s ^ ^2s-i,2s+2 
(5 = 1, • • • ,g - 1) and B2q-i^2q ^2,2q-i, wMch gives 

— (t — A)'^e~*^^. These permutations are independent 
of each other. Our numerical study shows that when 
/c = and 7r/g, the sum of all terms generated by 
the first kind of permutations is zero except the full- 
permutated ones (if there is) without any —jig left. So 
the first kind contributes a term of for odd g^'s or 
(A^ - fy/^ + (A^ - t2)9/2e-2^^9 for even g's. Thus we 
have Pf[5(0)] = - riLi - [(t + A)^ + (t - A)^] for odd 
^'s or nLi - [(t + A)^/2 - (-l)^/2(t- A)^/2]2 for even 
g's, and Vi[B{^/q)] = - [iLi + [(t + A)^ + (t- A)^] for 
odd ^'s or nLi + [{t + A)^/2 + (-l)^/2(t - A)^/2]2 for 
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FIG. 5: The spectral functions and the disorder effects. Here, 
as example, we choose a = 1/2, ^ = to plot (a) the 
momentum-resolved spectral function, A{k^uj)^ and (b) the 
corresponding total spectral function, A{uj)^ without disor- 
der, i.e. = 0. (c) shows that the total spectral function, 
A{uj)^ varies with the strength of disorder. The other param- 
eters are L — 250, t = A = F = 1 under OBC. We have done 
500 realizations to calculate the disorder effects. 



even g's. Therefore, we can get the general formula for 
the boundaries: IflLi Ms| = (^ + + (^ - for odd 
q\ and - [iLi = + ^Y + - ^Y + 2(A2 - t^y/^ 
(rc=i < 0) or nLi ^s = it^ AY + (t - A)^ - 2(A2 - 
^2jg/2^ (ns=i ^) even g's. The special cases 

of a = 1/2 and 1/3 can be derived from these general 
formula. 

Spectral function.- As the Majorana edge states corre- 
spond to the zero mode solution protected by the pres- 
ence of an energy gap, one would expect a zero-energy 
peak appearing in the corresponding spectral function. 
The momentum-resolved spectral function with momen- 
tum k and energy uj {h = 1) is defined as A{k,uj) = 
-^lmGr{k,uj), where Gr{k,uj) = J^^Gr{k,t)e'^^-^^^dt 

and Gr{k,t) = -iO{t) (G| {cfe(t),4(0)} \G) is the single 
particle (retarded) Green function with \G) being the 
ground state of the system, and Ck{t) = e^^^c/^e"*^^ is 
the fermionic annihilation operator in the moment space 
with Ck = -^^j Cje~'^^^ . And the corresponding to- 
tal spectral function is A{uj) = A{k^uj). As shown in 
Figin^a) and (b), an obvious zero-energy peak is observed 
for the system with OBC. On the contrary, we find no 
zero-energy peak in the system with PBC. 

The Majorana edge states are expected to be immune 
to local perturbations. To explore the effect of disor- 
der, here we add a random on-site potential, i^D = 
^iWiclci, to the Hamiltonian ([1]), which usually leads 
to Anderson localizations, where Wi is uniformly dis- 
tributed in the range of [—W/2,W/2] with W being the 
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disorder strength. Fig. [5jc) shows that even though the 
energy gaps between the zero mode and the excitation 
modes are smeared by disorder, the zero-energy peak is 
stih notable, although it is becoming shorter as well. It 
again indicates that the MFs appearing in our system are 
robust against local perturbations as expected. 

Summary.- In summary, we find that in the density 
modulated p-wave SC wires, by tuning the phase param- 
eter, the zero mode Major ana edge states emerge in some 
intervals, and some of them are strong enough, indepen- 
dent of the strength of the periodic density modulation. 
The appearance of the MFs demonstrate the Z2 topologi- 
cal nature. After calculating the Z2 topological invariant, 
we get the phase diagram for the transition from topolog- 
ically nontrivial phases to trivial phases. By this model, 
we supply a good platform to have possible schemes of 
searching for MFs. At last, we also give an evidence in 
the spectral function where a zero-energy peak appears 
under OBC, even subject to strong disorder. The spec- 
tral function is possible to be experimentally detected by 
the photoemission spectroscopy. 
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